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Abstract 
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^+ , By using the product ansatz as an approximation for the two-baryon system we investigate the 

■"".^ ■ isoscalar nucleon-nuclcon spin-orbit potential in an extended Skyrme model including both fourth- 

O ,' and sixth-order terms. As it is the case for the Skyrme model, we still obtain the wrong sign for this 

f^ i' interaction. Nevertheless, concerning the order of magnitude, the extended Skyrme model provides 

r^ , a better agreement with phenomenological potentials, as compared to the standard one. 
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The Skyrme model ^, recognized as the simplest chiral realization of QCD at low energy and large 
Nc Q, has been extensively used to describe static and dynamical properties of baryons. However, it 
has been shown recently |y, Q, through the study of some properties of the nucleon, that one should 
not restrict oneself to the standard Skyrme model for the description of low-energy hadron physics, 
but consider extensions of this model including higher order terms in powers of the derivatives of the 
pion field. In this letter we wish to study the effects of such additional terms on the dynamical part 
of the nucleon-nucleon interaction. We examine the influence of the sixth-order term generated by uj- 
meson exchange [^ on the nucleon-nucleon spin-orbit force, where we focus on the isospin independent 
part of the interaction. Several attempts to extract the isoscalar N-N spin-orbit interaction from the 
standard Skyrme model have given the wrong sign M, M. Riska and Schwesinger P] claimed that the 
inclusion of the sixth-order term leads to the right sign. However, these authors used a value for the 
corresponding parameter which is in conflict with the experimental uj ^ tt'j width. 

1 . The effective Lagrangian density corresponding to the extended Skyrme model we use can 
be expressed in terms of an SU(2) matrix field U which characterizes the pion field. It reads 

£ = ^ Tr {d^Ud^U+) + :^ Tr { [{d,U)U+ , {d,U)U+]' } 



1 B^ P 

2m2, ^ 4 



--^B^B^' + ^m^ Tr {U + U+ - 2) 
2 mi 

where B^ is the baryon current 



^^ = ^^'''"''T^ i {d,u)u+{d^u)u+{dpU)u+ }. (2) 

The first term in Eq. (||) corresponds to the nonlinear cr-model, /tt being the pion decay constant. The 
second term, which is of fourth-order in powers of derivatives of the pion field and parametrized by 
the dimensionless factor e, was introduced by Skyrme in order to stabilize the soliton. The third term 
is of order six in the derivatives and can be derived from a local approximation of an effective model 
with uj mesons Q. The constant rrii^ is the cj-meson mass (782 MeV) while j3i^ is a dimensionless 
parameter related to the uj ^ t^^ width. The last term in Eq. (^ which is proportional to the square 
of the pion mass m-„: (139 MeV) implements a small explicit breaking of chiral symmetry. 

For the system of two interacting solitons we use the product ansatz as suggested by Skyrme @ . 
We also introduce rotational dynamics to obtain the appropriate spin and isospin structure |10|. Thus, 



the field configuration of the two-nucleon system separated by a vector r reads: 

UB=2^U2{^,r,A,B)=AUH{^-r/2)CUH{^ + r/2)B+ , (3) 

where C = A^B and A and B are SU(2) matrices. To carry out a simultaneous quantization of the 
relative motion of the two nucleons and the rotational motion we need to treat r, A and B as collective 
coordinates. Hence we make all these parameters {y,A,B) time dependent. In Eq. (|^, Uh is the 
commonly used SU(2) matrix for a single soliton with the hedgehog ansatz: 

C/h(x) = exp[ if.xF(|x|) ] , (4) 

where i^(|x|) obeys the usual boundary conditions for winding number one, and the Tq's are the Pauli 
matrices. The notation x means x/|x| . We have to note that the product configuration @ is an 



approximation which has inconveniences and advantages |11]. We choose it for its relative simplicity 
as compared to other two-baryon field configurations which can be found in the literature. The region 
of validity of the product ansatz corresponds to a separation r much larger than 1 fm. 

2. The spin-orbit potential will emerge due to a coupling between the relative motion and the 
spins of the two nucleons. So we have to calculate the kinetic energy 

K = J d^x ICiUi) , (5) 

where /C(f/2) is that part of the Lagrangian (Q) involving only quadratic time derivatives, in which U 
is given by Eq. ^. The full expression of Eq. (|5|) with respect to {r,A,B) and their derivatives is 
very complicated and contains several terms which do not need to be written in this letter. They will 



be reported on a forthcoming publication |12]. The isospin independent part Kq of the kinetic energy 
(H) reads 

Ko = ^QMoQ , (6) 

where 

Q = (r,w+,w_). (7) 

In Eq. (0), w-t are the sum and the difference of the rotational velocities respectively: 

i 



w± = -- (Tt{tA^A) ± Tr(f5+S)) , 



and Mo is the 9x9 mass matrix 

/a{r)I -I3{r) fl \ 

Mo = /3(r) 4 7(r) / (8) 

\ 7(r)I/ 

where / is the 3x3 identity matrix and ^ij = eij/.i'k- The functions a,/3 and 7 depend only on the 
relative distance r = |r| and are determined from the chiral function F (see Eq. (B)). We have 
to mention that the mass matrix Mo contains other terms, such as a coupling between rotational 
velocities or a radial motion term, which contribute in principle to Kq. However, we find these terms 
very small and thus neglect them |0, ^. In this sense, the matrix Mo given by Eq. (|8|) should be 
considered as the dominant contribution to the isospin independent kinetic energy (^). 

Before identifying and extracting the spin-orbit potential from Eq. (^ one has to treat carefully 
the conversion from velocities to canonical momenta. Indeed, one has to invert properly the mass 
matrix Mo in order to move from a Lagrangian formalism to a Hamiltonian one. In terms of the 
canonical momenta 

P = MoQ = (Pr,s,s_) , 

where Pr is the conjugate momentum to r, and s = 1 , s_ = the conjugate momenta 

to vif+, w_ respectively [^, |l^ ((?i/2 being the spin of the nucleon i, i = 1,2) , the kinetic energy (^) 
becomes: 



Ko = ^PMo'ip . 
It is straightforward to invert the mass matrix (pi). It's expression reads: 

f ii H \ 
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where ^ has been defined above and fij = rifj. By developing Eq. (|9|) we obtain 



2i^n 



1 



a7 — /3^ 



7 P: 
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(11) 



where I = r x pr is the angular momentum. Now we can easily extract from Eq. ( pj]) the isoscalar 
spin-orbit force with the result: 



Vis{r) 



+ 



r (a7 — /32) 



(12) 



3. In order to compute the spin-orbit potential (|T^ we have to take a set of parameters of 
the model, namely, fj^, e and /J^;. The physically more sensible manner to choose the parameters is to 
remember that the Skyrme model involves only meson fields, and where baryons emerge as topological 
solitons. Therefore, since it is a meson theory, the parameters of the model have to be fixed by fitting 
to the low-energy meson observables and not to baryon data as it has been done in Ref. |^. The 
experiment yields /t^ = 93 MeV for the pion decay constant. The dimensionless parameter e can 
be determined from chiral perturbation theory |14]. The best set of chiral low-energy constants has 



been determined by Riggenbach et al |l^], in analyzing the K^ decays. As a consequence, one finds 
e = 7.1 ± 1.2 Q. Finally, the last parameter Pui is obtained by fitting to the w — > 7r7 width, yielding 
Pui = 9.3 flH]. We have to mention that this set of parameters yields values of the nucleon mass 
(including the Casimir effect), the A-A^ mass splitting, the axial- vector coupling constant, the Roper 
resonance, etc., close to the experimental situation ||3|, |^. 

Our first result is that the isoscalar spin-orbit force obtained with these parameters is repulsive and 
thus has the wrong sign as compared to the experimental situation. This disappointing prediction of 
Skyrme models is still not explained until now. In order to compare our results to the data we consider 
therefore —Vis instead of Vis (cf., Eq. (|12|)). Of course one can force the interaction to be attractive 
by increasing the value of the parameter P^ and thus the sixth-order term dominates the second- and 
fourth-order terms. This was done by the authors of Ref. |8[ who used a value of Pui = 22.4. However, 
this way to proceed is not acceptable since this value of P^ disagrees stongly with the experimental 
(J — > 7r7 width. 

We plot in Fig. 1 the behaviour of —Vis (0) with respect to the nucleon- nucleon separation r for 
[f-K = 93 MeV, e = 7.1,(3^ = 9.3] and compare it to the standard Skyrme model prediction ((3^^ = 0). 
We also plot in the same figure the corresponding part of the phenomenological Bonn potential [p^ . 
From Fig. 1 we obviously see that the inclusion of the sixth-order term in the Lagrangian (|l|) improves 
considerably the prediction of the isoscalar spin-orbit force in the region of validity of the model ,i.e., 
r > 1 fm. It is expected that the agreement can not be perfect due to several reasons. One of them is 
the use of the approximation of the product ansatz (y) to describe two-interacting nucleons. A second 
reason is that the model used here (IT]) should be regarded as a minimal extension of the Skyrme 
model. Concerning the region of smaller r (not shown on the figure) the prediction of the model 
differs completely from the phenomenological potential because this region corresponds to processes 



of large momentum transfers where one should not use effective Lagrangians but perturbative QCD. 

4. In this letter we have derived the isospin independent spin-orbit interaction from the 
extended Skyrme model including, in addition to the nonlinear a model, fourth- and sixth-order 
terms in powers of the derivatives of the pion field. Concerning the order of magnitude, we find 
a considerable improvement as compared to the Skyrme model prediction and thus, confirm the 
conclusions of Refs. |^, ^ . Indeed, the standard Skyrme model is insufficient for an accurate description 
of baryon phenomenology and one should generalize this model by including higher order terms in the 
chiral Lagrangian. In spite of that, our major result is that, even with the sixth-order term present, we 
still get the wrong sign for this force. One should try to understand why the standard Skyrme model as 
well as the extended one (|^) give the wrong sign and not mask the problem by considering irrealistic 
values of the parameters |^. This difficulty with the isoscalar spin-orbit force seems to be similar 
to that of missing central attraction. The problem concerning the lack of attraction in the central 
potential has been solved by considering effective Lagrangians which incorporate low mass mesons 
with finite mass (the scalar meson having an important role) |jl^. Maybe one has to investigate in 
this way in order to correct the anomaly of that sign [^ . 
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Figure captions 



FIG. 1. The isoscalar spin-orbit potential (12) with a minus sign in the case where the parameters 
fit to low energy meson data. The dotted line corresponds to the Skyrme model (i.e., f3uj = 0), while 
the full one corresponds to the extended model with P^^ = 9.3. The dashed line is the corresponding 



term in the Bonn potential |17]. 
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